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SUMMARY 

Equations relating the critical temperature and ion current density 
for surface ionization of cesium on tungsten are derived for the cases 
of zero and finite electric fields at the ion-emitting surface. These 
equations acre used to obtain a series of graphs that can be used to 
solve many problems relating to ion-rocket theoretical performance. The 
effect of operation at less than space -charge -limited current density 
and the effect of nonuniform propellant flux onto the ion-emitting sur- 
face are also treated. 


INTRODUCTION 

The potential of the ion rocket as a space propulsion engine (for 
example, refs. 1 and 2) has led to a sizable research and development 
program on this engine. This work is being conducted in both the NASA 
Lewis Research Center (ref. 3) and in several laboratories in industry. 

A method for estimating the performance to be expected with various ion- 
rocket configurations at different operating conditions is therefore of 
importance. 

Figure 1 depicts one simple ion -rocket configuration wherein a con- 
tact ion source is employed. Many of the concepts and terminology that 
follow can be understood by reference to this figure. 

The type of ion source receiving the greatest research attention 
for possible ion-rocket application is the contact ionization of alkali 
metals, primarily cesium, on heated surfaces having a high electron work 
function, primarily tungsten. This report presents a graphical method 
for estimating the performance of those ion rockets that employ contact 
ionization of cesium on tungsten surfaces. The relations on which the 
graphs are based were obtained from theory and from analysis of reported 
data on contact ionization. The various assumptions and limitations 
that apply to the method are listed in a subsequent section of the report. 
The report also considers the effect on efficiency of nonuniform current 
density and of operation at current densities below the space-charge- 
limited condition. 
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SYMBOLS 

ion-emitting area, 
electric-field strength, v/m 

distance of propellant supply tube to emitter surface, m 

specific impulse, sec 

ion current density, amp/m^ 

average ion current density over total emitter area, amp/m^ 
ion current density for ideal diode, *tmp/m^ 
maximum local ion current density, amp/m^ 

Boltzmann constant 
accelerator length, m 
ion mass, kg 

power input to accelerator electrode, w 
ion beam power, w 

power required to emit electrons, w 

power input to ion rocket from general ,or, w 

power required to vaporize and ionize propellant, w 

summation of power losses from ion rocket, w 

power radiated from ion emitter, w 

energy loss by radiation from ion source, ev/ion 

electronic charge, coulomb 

ion emitter temperature, °K 

voltage , v 

voltage at ion emitter, v 
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AVa voltage drop across accelerating gap, v 

AVi p overall voltage drop, v 

W e width of emitter surface, m 

X fractional distance from emitter edge 

Z ratio of beam power to radiated power at two conditions 

(defined by eg. (3l)) 

e thermal emissivity 

Gq permittivity of free space 

T] ion-rocket energy efficiency 

r] s energy efficiency of ion source (uniform flux) 

Tj ^ energy efficiency of ion source (nonuniform flux) 

cp potential energy, ev 


EQUATIONS FOR ION-ROCKET EFFICIENCY 
Ion-rocket energy efficiency is given by; 


n = 




P B 


loss 


r B 


gen 


a) 


The various components of the term ^^loss ^ or 311 ^ on roc ket employ- 
ing a contact ion source are as follows; 


Pr The radiated power from the hot tungsten surfaces. 

This is by far the largest part of ^Pq oss in most engines using con- 
tact ion sources. For a given P^, the value of P^ attainable is 

dependent on the uniformity of the propellant flux onto the ionizing 
surface; this will be further discussed. 

Pp ro p The power required to vaporize and ionize the propel- 
lant. This is usually negligible: ?p ro p < 6J, watts, where J is 

total current in amperes. 
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Select The P ower required to emit tie electrons for beam 
neutralization. The value of this term can range from very small to 
very large values relative to the other power losses, depending on the 
electron emitter configuration and work function (hence, temperature). 

The power required to maintain electric potential on 

the engine electrodes. This term must be very small in any engine having 
reasonable operating life, since ion impingement on the electrodes or 
electron emission from the electrodes would result in '•severe erosion 
and/or overheating of engine parts. 

Of the various component terms in ^^Pq oss , only P R is included 
in the computations and graphs that follow. The term P R is associated 
with the ion source, while the other ^2 P 1 os 3 component terms are as- 

sociated with other parts of the engine. In the strict sense, therefore, 
the following discussion deals with the ion-source energy efficiency q s : 


% 


P B + P R 


( 2 ) 


As just noted, the other terms in ^P loss nay prove to be negligible 
for many ion rockets, in which case 

T1 S ~ H ( 3 ) 


Many of the equations that follow apply onlj to q s $ others apply also 
to q, and in such cases the appropriate syirbols are used. 


t? 


GRAPHS FOR I ON -ROCKET PERFORMANCE 

Most of the equations for ion-rocket theoretical performance are 
well documented (e.g. , ref. 3); derivations of these equations are not 
repeated herein. Except when otherwise noted, the equations that fol- 
low are included in, or can be easily derived from, those of reference 3. 


Performance with Uniform Prcpellant Flux 

For cases where propellant flux onto the ionizing surface is uni- 
form, figure 2 permits a graphical de terminal ion of various performance 
parameters, including q s . Figure 2 contains four page-size graphs 

that can be used to determine ion-rocket theoretical performance char- 
acteristics. Using figure 2(a), the overall voltage AV T required to 
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yield a given specific impulse is determined. Using figure 2(b), for a 
given accelerator length L and accelerating voltage AV^, the ideal 

current density Jideal can rea( ^ from the abscissa. This value of 
j ideal can used to enter figure 2(c), where j and emitter surface 
temperature T serve to determine Q the radiant-energy loss from 

the emitter in units of electron volts per ion. The lower boundary 
curve results from the minimum permissible temperature of operation for 
the emitter surface governed by the current-critical temperature rela- 
tion, equation (7). The value of can then be used to enter fig- 

ure 2(d), where Q R , together with the overall voltage drop AV T pre- 
viously obtained, serves to determine r] s j rj s is the ion-source energy 
efficiency for the ion rocket giving the I g specified in figure 2(a). 

A graph containing all four parts of figure 2 on a single sheet is in- 
cluded in addition to the page-size graphs of parts (a) to (d). This 
arrangement of the four parts facilitates their use. 

The various parts of figure 2 are subject to the following assump- 
tions and limitations: 

(1) The propellant is cesium. (Figs. 2(a) to (c). ) 

(2) The ion-emitting surface is tungsten with emissivity given by 
figure 3. (Fig. 2(c).) 

(3) Operation is at the critical temperature defined by equation 
(7). (Fig. 2(c). ) 

(4) Current density is equal to that for space-charge-limited 
current from an ideal plane -parallel diode of spacing L. In some cases 
this assumption can be modified by techniques to be discussed in a 
later section. (Fig. 2(b). ) 

(5) Thermal radiation occurs only from the ion-emitting surface. 
This implies complete radiation shielding on the upstream face of the 
emitter and zero shielding downstream. (Fig. 2(c). ) 

(6) Space charge limitations occur only at emitter. (Figs. 2(b) 
and ( c ) . ) 

The essential equations used in obtaining the graphs of figure 2 
are as follows; 

Graph ( a ) : 


I s = 123 VAV t 


(4) 
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Graph (h): 


j = 0. 475X10" 8 


(AV a ) 3 / 2 


(5) 


Graph ( c ) : 

( 1 ) % - 10 - 28 J (so ) 4 (6) 

(2) Emissivity of tungsten as a function of temperature as given 
by figure 3 


( 3 ) log 10 j = 11. 3 


13,200 

T 


(7) 


(Derivation of this equation discussed in later section of report) 


Graph ( d ) : 


T) = 


fB_ 

P B + P 


R ^R ^R 

+ % 1 + z% 


(s) 


If aperture and space charge effects are known for a particular 
ion rocket so that j / j ideal' this correctiorL can 136 applied on the 
abscissa of figure 2(c) before determining Q, R . For grid electrodes 
the values of j/j ideal can es ' fcijnated hy the method of reference 4 
(see also ref. 5). For other electrode configurations, j/j ideal can 
often be determined from limited test data on the particular engine 
under consideration. 


Performance with Nonuniform Propellant Flux 


In some ion rockets a nonuniform propellant flux onto the ion- 
emitting surface may be encountered. With porous tungsten emitters 
this might result from nonuniform porosity. In other contact ion 
sources it might result from the pattern of flu: out of the propellant 
vapor supply system. One example of nonuniform propellant flux is shown 
in figure 4$ this is the computed flux pattern :‘or the engine depicted 
in figure 1. 

Where a nonuniform propellant flux exists, the ion emitter tem- 
perature must be maintained at the critical temperature corresponding 
to the local j max - Otherwise, an adsorbed ces .urn layer will build up 
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at that location with the resulting emission of a large percent of 
neutral cesium atoms. Thus, while T and. P R are fixed by J max , 

the value of Pg is reduced from the value that would be obtained if 
J we re uniform across the emitter with j = j max . This reduction in 
Pg results in a lower value of rj s and tj. 


In figure 4 the area under the curve is proportional to j av , while 
the total rectangular area is proportional to J max - Now, if tj repre- 
sents the energy efficiency that would be obtained with uniform propel- 
lant flux and T]g represents the efficiency with nonuniform flux, then 


1 = 


J max 


^max ^ ^ 

A AV + Jjloss 


( 9 ) 


and 


^ = 


Jav A AV 


J av 


AAV + 


E 1 


loss 


(io) 


Now, values of rj (actually t] s ) are obtained from figure 2, and 
simple correction factors are desired to use with these values to ob- 
tain T]g and T] g ^g for cases where propellant flux is not uniform: 





( 11 ) 


and 


Hs.S 

^s 


A plot of rig/ r) and ti s g/rjg 

rocket would have efficiency t| or t] s with uniform propellant flux 

(fig. 2), the efficiency of that same ion rocket with nonuniform flux 
characterized by some value of j av /j max is given by appropriate 

values of rig/ri or t] s /t} s from figure 5. 


Jav/ J 


a v/ d max 


&-■) 


1 J 7 ls + 1 

"max 

is included in figure 5. 


( 12 ) 


If a given ion 


A comparison of the efficiency for uniform and nonuniform propel- 
lant flux can also be made readily by use of figures 2(c) and (d). Th 
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efficiency r\ is read for j max as described previously. The effi- 
ciency Ti s is obtained as follows: From the point in figure 2(c) cor- 

responding to j max and critical temperature (lower line), go up a 
constant -temperature line to j av - This poinn gives the appropriate 
value of Qr, and the efficiency tj s is then read from figure 2(d) at 
this value of Q, R and the appropriate AV T . 

CRITICAL TEMPERATURE - CURRENT RELATIONS 
IN CESIUM-TUNGSTEN SYSTEM 
Zero Electric Field at Emitter 

A different current - critical temperature relation (eq. ( 7 )) is 
used in this report from that used previously (refs. 3 and 5). This 
equation was derived from the data of table I of reference 6 and is 
recommended as the relation to use for space -charge -limited flow con- 
ditions, that is, for zero electric field at the ion emitter surface. 

The equation relating current and temperature that previously was used was 
obtained from data at other than space-charge-limited flow conditions. 

A comparison of the equation used herein with the earlier equation at 
current-dens ity levels of interest for ion-rccket engines indicates a 
temperature difference of around 150° K between the two equations. This 
difference may make considerable difference in the power radiated by 
the ion emitter and hence in the calculated ion-rocket efficiency. 

The zero-field relation presented hereir was obtained by plotting 
the total evaporation rate (sum of the ion axd atom currents) against 
the fraction of a monolayer adsorbed on the surface for several constant 
surface temperature values (fig. 6). A discission of the behavior of 
the cesium-tungsten system as governed by these total-evaporation-rate 
curves may be found in reference 7. In particular, it may be noted that 
a plot of the current at the minimum in these curves against the tem- 
perature will yield a current - critical temperature relation. The 
equation so obtained is 


i°gio J 


11.3 


15,20 0 

T 


(?) 


A plot of this relation is shown in figure 7. 
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Electric Field at Emitter 

The existence of an electric field at the ion-emitting surface 
will affect both the current -voltage relation and the current - critical 
temperature relation. The question arises as to whether there is an 
optimum condition for operation, that is, one having the lowest ratio 
of radiated power P^ to ion beam power Pg. The discussion in this 

section is directed toward an answer to this question. 

The equation for current flow between parallel planes when an 
electric field E exists at the emitter plane can be expressed 

L = -^[(bVAV - 2E 2 ) (b-y/AV + E 2 ) 1 / 2 + ZE^] (13) 


where 



A plot of equation (13) for cesium in the range of AV = 40,000 
volts and j = 200 amperes per square meter is shown in figure 8. 

The effect of the existence of an electric field at the ion- 
emitting surface will be to reduce the effective potential barrier for 
the escape of ions and hence to increase the ion current flow, similar 
to the Schottky effect for electrons (ref. 8), illustrated pictorially: 



The symbol cpQ represents the potential energy barrier to the 

escape of an ion in mutual interaction with the electrical image force. 
The amount of the reduction in the potential barrier due to the pres- 
ence of the field is given by the relation 
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( 14 ) 


where Ztfp is expressed in electron volts and E in volts per meter. 
Simplifying, this expression becomes 

Ap = -3. 8X10" 5 E 1 / 2 (15) 

In order to illustrate the effect of the field on ion current flow, 
equation (7) can be written in the form 

j = 10 11 ' 3 e- 2 -62/K (16) 


In this form, from the exponential Boltzmann factor, the potential 
energy barrier to escape of ions at zero field would be interpreted as 
2. 62 electron volts. 


The change in current flow caused by the field may then be given 




jE,T 0 - jE=0,T 0 


-Ap/kT 


(17) 


where 


j e=0,t 0 

JE,T 0 
and Z^p, 


current at zero field 

current at field E and same temperature 
the change in potential barrier, is given by equation (15). 


An equation for current - critical temperature with other than 
zero electric field can be generated by using equations (15) and (17) 
in conjunction with the zero-field ion current relation. 


Incorporating equations (15) and (17) i ito (16) yields the fol- 
lowing relation for current - critical temperature at other than zero 
field at the emitter: 


j 



(18) 


. 11. 3 - S #< 1 - 1 - 45X10 ‘ 5 V®) 

J = 10 e 


(19) 


or, in logarithmic form, comparable with equation (7), 

log 10 j = 11.3 - + 0.192 ^ (20) 
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This relation, encompassing the effect of electric field on the current 
critical temperature relation, may he compared with the experimental 
data of reference 9. In reference 9 data are presented on the effect 
of electric fields up to 2X10® volts per meter on the critical temper- 
ature for constant cu rrent. Valu es of dT/d-^/E ranging from about 
-0. 01° to -0. 015° K/-\/ volt/me ter were obtained. 

By differentiating equation (20) at constant current, the following 
expression for (dT/d-/E)j is obtained: 

/ dT \ = _ 1. 45X10 ~ 5 T 

1 - 1. 45X10" 5 -y/E 

If E is limited to about 10^ volts per meter, this relation is 
approximately expressed 


dT r 

rz -1. 45X10 ’ d T (22) 

a-/E 

The range of temperatures investigated in reference 9 was about 700° 
to 900° K. In this temperature range, values of equation (22) vary 

from about -0. 01° to -0. 013° K/ -/volt/meter, which is in agreement with 
the experimentally determined values. 


Efficiency with Other than Zero Field at Emitter 


To determine whether there is an optimum condition of operation, 
that is, one in which the ratio of "beam power P B to radiated power 
% is a maximum, the ratio P b/ P R at off -space-charge conditions can 
be compared with that at the space -charge-limited condition 

p b/ p r _ p b/ p b , 0 
P B,o/ P R,0 P r/ P R,0 


J/J 0 


(^AqT q)- 


= z 


(23) 


where the subscript 0 represents the space-charge-limited condition. 

The actual ratio of efficiencies tj/tJq will not vary as much as 
the function in equation (23), since 


JL = Z 

T} 0 1 + T] 0 (Z - 1) 


(24) 


With tungsten, in the temperature range of interest for surface 
ionization, the emissivity is proportional to temperature, so that 
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Z = 


j/jp 

(T/T 0 ) 5 


Equation (25) can "be rewritten: 



(25) 


(26) 


The ratio Z may he readily evaluated for a constant accelerator 
spacing L and constant accelerating voltage AV. As illustrated in 
the following sketch, is the space -charge -limited current flow. 



When the current is lowered hy an amount Aj_. a field E is present 
at the emitter (determined from eq. (13) or fig. 8). The required 
critical surface temperature is lowered, both because of the field and 
because of the lowering of the current: 

AT = AT B + «Pj - - (g)Aj (27) 

Substituting from (22) for 5 t/B-\/E and froai equation (16) (by first 
differentiating) for ST/Sj gives 

^ = -1. 45X10 -5 VE + — n ^sl (28) 

T V 2. 52 j 

and thus equation (26) can be evaluated by the relation 


+ (Aj/jp) 

[l - 1. 45X10-5 yE + ^ (^-)] 5 
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Values of Z have been evaluated at several sets of conditions 
varying from AV = 10,000 to 50,000 volts and currents from 200 to 1000 
amperes per square meter and are shown in figure 9. According to these 
curves, more efficient operation could be attained at current densities 
a few percent below space-charge-limited values. At high current den- 
sities the region of greater efficiency for off-space-charge operation 
is greater. 

It should be noted, however, that (l) the percentage improvement 
in efficiency for off-space-charge operation is small, even less than 
the values of Z (note eq. (24)); (2) the surface temperature lowering 
permitted to take advantage of this off -space -charge operation advantage 
is small - of the order of 50° C or less. Such close control of emitter 
surface temperature may be difficult. 


CONCLUDING REMARKS 

The graphs of ion-rocket performance parameters presented herein 
constitute a means for rapid determination of theoretical engine per- 
formance characteristics. The equation relating emitter temperature to 
space -charge -limited current density, on which the performance graphs 
are based, was derived from experimental data obtained with pure tung- 
sten ion emitters in small-scale apparatus. NASA data on an experi- 
mental ion rocket (ref. 10) indicate that higher emitter temperatures 
are required in that engine than those indicated by this equation. 

This discrepancy may be due to the presence of surface impurities on 
the engine ion emitter. Such impurities, if present, might also alter 
the thermal emissivity of the emitter from that of pure tungsten. Thus, 
in ion-rocket experiments, the presence of impurities on the ion-emitting 
surface can result in radiation power losses that are greater than those 
computed from the equations given in this report. However, It should 
be possible to attain the engine performance indicated herein through 
the use of inerting and bake-out procedures to keep the ion emitter 
free of contaminants. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, May 12, 1960 
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(a) Specific impulse. I g = 1‘ 3-y/AVrp. 

Figure 2. - Graphs for determining ion-ro< ket -engine performance. 
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Figure 2. - Continued. Graphs for determining ion-rocket-engine performance. 
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Fractional distance from emitter edge, X/W e 


Figure 4. - Typical nonuniform current distribution pattern for 

engine shown in figure 1* Distance of propellant supply above 
emitter to emitter width, h/W e , 0.3. 
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Figure 5. - Effect of nonuniform propellant flow on ion-rocket 
energy efficiency. 
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Fraction of monolayer adsorbed on surface 


Figure 6. - Total (ion + atom) evaporation of cesium on tungsten for zero 
electric field at surface. 
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Figure 7. - Current - critical temperature relation (eq. (?)) for cesium on tungsten. 
log 10 j = 11.3 - 13,200/T. 
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Figure 8. - Typical plot of voltage gradient at emitter with 
changes in accelerator length, voltage, and current (eq. (13)). 
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